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^ ' Abstract. A method is presented for obtaining rigorous error estimates for approx- 

, O , imate solutions of the Riccati equation, with real or complex potentials. Our main 

j /' j ■ tool is to derive invariant region estimates for complex solutions of the Riccati equa- 

tion. We explain the general strategy for applying these estimates and illustrate the 
method in typical examples, where the approximate solutions are obtained by glue- 
ing together WKB and Airy solutions of corresponding one-dimensional Schrodinger 
equations. Our method is motivated by and has applications to the analysis of linear 
i-i^ . wave equations in the geometry of a rotating black hole. 

^ ■ 1. Introduction 

-)— » 

The time-independent one-dimensional Schrodinger equation plays an important 
role in both mathematics and physics. Except in special situations, it is impossible 
to find exact solutions, and thus one resorts to approximation techniques. The most 
Cn ' common method is to take the WKB wave functions away from the zeros of the poten- 

tial, and to glue them together with Airy functions, which approximate the solutions 
Q I near the zeros (see for example the textbook [9l Section 2.4] or the book [5j). Usually, 

'nJ" ■ the potential in the Schrodinger equation is real-valued. However, there are situations 

(like systems with dissipation) where the Schrodinger potential is complex. 






^ 



^ 



C"^ I The shortcoming of most approximation procedures is that they do not give error 

^Q ■ estimates. To quote Froman and Froman [51 page 3] (1965), 

^^ . "In spite of the abundant literature there is still lacking a convenient 

method for obtaining definite limits of error. . . The problem of obtain- 
ing such limits of error is not only of academic or mathematical interest 
but is especially important in some physical applications. . . " 
5^ I Unfortunately, such a convenient method for obtaining definite limits of error is still 

lacking today. Indeed, to our knowledge, despite the vast literature on WKB methods 
(see [1] and the references therein), only a few papers are concerned with rigorous 
error estimates. As typical examples see [71|8], where L^-estimates are derived in the 
time-dependent setting for short times, as well as [6|, where estimates are obtained 
in the setting of multiple tunneling. Here we shall develop a general method for 
obtaining rigorous pointwise error estimates for any given approximate solution in the 
time-independent setting, even in the case when the Schrodinger potential is complex. 
We were led to the study of error estimates for solutions of one-dimensional Schrodin- 
ger equations from problems arising in general relativity. In particular, when analyzing 
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the scalar wave equation in the Kerr geometry, Carter's remarkable separation of 
variables [1] gives rise to a radial and angular ODE. These equations can both be 
transformed into the Schrodinger form, which are coupled through two parameters u 
and A. For estimating the global behavior of the solutions uniformly in uj and A, we 
were led to the idea of using invariant disk estimates [3, 2j. When attempting to extend 
our analysis to gravitational waves in the framework of the Teukolsky equation |10j . 
one of the difficulties that arises is that the potentials in the separated equations turn 
out to be complex. In the present paper, we obtain invariant disk estimates in a 
general setting, which are applicable to one-dimensional Schrodinger equations with a 
real or complex potential. The previous estimates in [3l [2] will arise as special cases. 
In particular, our method gives rigorous pointwise error estimates for the standard 
approximations by WKB and Airy functions. But the method works just as well for 
any other approximate solutions. The aim of the present paper is to introduce the 
method from a general point of view, and to illustrate it in typical examples. For an 
application to eigenvalue and gap estimates as well as to problems involving parameters 
we refer to [3] and [2]. 

We now set up the framework and fix our notation. The one-dimensional Schrodinger 
equation reads 

' 6? 



..= +''(-' 



4){x) = \4>{x) 



Absorbing the parameter A into the potential, we rewrite the Schrodinger equation as 

^cA = V{x)(p{x), (1.1) 

where the new potential V = V — \ may depend on additional parameters. This is 
a linear second order ODE. Thus the general solution (j) can be written as a linear 
combination of two linearly independent fundamental solutions 0i and 02, 

(j) = a(Pi+b(j)2- (1.2) 

By choosing the coefficients a and b appropriately, one can then satisfy suitable bound- 
ary conditions, if so desired. Setting 

y = ^ , 1.3 

the function y satisfies the Riccati equation 

y' = y-y' . (1.4) 

Our strategy is to estimate y in the complex plane. This requires a brief explanation, 
because in the case when V is real it is natural to also choose 4>, and thus y, to 
be real. This real function <j) will in general have zeros, and y will be singular at 
each of these zeros. These divergences make it difficult to estimate y. We avoid 
this problem by considering a linear combination (p = (pi + i(j)2 with real fundamental 
solutions (pi and (/)2. Then, since the zeros of the fundamental solutions never coincide, 
the function (j) never vanishes, and thus the corresponding function y is globally well- 
defined. 

In the more general case when the potential V is complex, the situation is more 
difficult because complex solutions of the Riccati equation ()1.4p may "blow up in finite 
time". Nevertheless, we can hope to avoid singularities of the Riccati equation by 
choosing the initial condition 2/(0) suitably. 
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In this paper, we focus on the analysis of solutions of the Riccati equation for 
specific initial conditions y(0). This is justified because our estimates also give rise 
to corresponding estimates for general approximate Schrodinger wave functions, using 
the following procedure. Suppose that a specific y is known, up to a small error term. 
Then integrating (jl.3p . one obtains one fundamental solution (j)i{x) = exp(/ y) of 
the Schrodinger equation. The other fundamental solution can readily be obtained 
by integrating the Wronskian equation (l)'i(j)2 — (j)i4>'2 = const. We thus obtain the 
general solution ()1.2p . Keeping track of the error in this procedure, one can derive 
error estimates for any Schrodinger wave function. 

2. The Flow of Circles for a Constant Potential 

In this section we consider the special case of a constant potential. This case can 
easily be analyzed in closed form, and the resulting flow in the complex plane will be 
helpful for the geometric understanding of the invariant disk estimates. Thus assume 
that F is a non-zero complex constant. We set ( = VV. Then the general solution of 
the Schrodinger equation (jl.ip is 

4){x) = ae^"" + he-^"" , x G M 

with a,b G C The corresponding solution of the Riccati equation is 

ae^^ — be~^^ (a — b) cosh(Ca;) + (a + 6) sinh(Cx) 



c 



ae"^^' + 6e~f^ (a + b) cosh(Cx) + (a — 6) sinh(Cx) 

(a - 6) + t{x) {a + b) 



(a + 6) + t(x) (a - 6) ' 
where t(x) := tanh(Cx). Setting yo = y(0), 

a 



2/0 = C 



a + b' 
we can express y in terms of yo, 

-rix) yo + C 
Note that the mapping yo ^ y is a fractional linear transformation. It has the fixed 
points y = ib(^ corresponding to the two solutions (j) = e^^^, respectively. Furthermore, 
we see that circles are mapped onto circles or straight lines. Thus by varying x, we 
obtain a flow of circles in the complex plane. To analyze this flow, let us assume that yo 
lies on a circle of radius Rq around ttiq, i.e. 

yo = "T-o + z with z = Rq e^'^ . 

Then for any fixed x, 

z + mo + rC 
Tz + TrriQ + C 
Thus for any m E C, we can write 

(C — T'm)Roe^'^ + (mo + r^^ — rmmo — m( c / e*"^ + 7 



m 



where the introduced quantities 

{C, — Tm)RQ (mo + rC^ — Tm,m,o — m,C, tRq 

rmo + C ' (C - Tm)Ro ' rmo + C 
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are independent of 99. For m to be the center of the new circle, the resulting expression 
for |y — ?7i| must be independent of (p. This gives the condition 7 = 5. We thus conclude 
that 

\y — m\ = R , 

where the center and radius are given by 



_ . {mo + TQiTmo + O-RlT Ro \t^ - 1| |C|' 

^ I ^19 n9 1 19 ' '^ 



\Tmo + C\^ - RI\t\^ ' \rmo + C\^-Rl\T\^' 

The denominator can have zeros, in which case the circle degenerates to a straight 
line. Whenever the denominator has no zeros, we can take the limit 2; ^ 00 to find 
that in the case ReC > 0, the circle asymptotically approaches the fixed point y = C> 
unless if we start at tjiq = —( with Rq = 0. Thus we can say that y = C is the stable 
and y = —( the unstable fixed point. Conversely, in the case ReC < 0, the fixed 
point y = C is unstable, whereas y = —^ is stable. The flow of circles is illustrated in 
Figured] (left). 

In the remaining case Re C = 0, we clarify the situation by computing the stationary 
circles. In this case, C and r are purely imaginary, and thus the above formula for R 
becomes 

R = Ro C 



r2-l 



^{Rl-\mo\^)-2TCRe{mo)-e- 
For a stationary circle, R must equal Rq, for all imaginary r (because the range of 
the function tanh((^x) is iM). Multiplying by the denominator, we get a quadratic 
polynomial in r, and comparing the coefficients, we obtain the conditions 

Remo = and |moP — -Rq ~ ICP • 
We thus obtain for any Rq >0 the two circles with centers 



mo = ±i^R'o + \C\^, (2.1) 

and an easy computation shows that these circles are indeed stationary under the flow. 
Since these circles enclose the fixed points (see Figured] (right)), both fixed points are 
centers. 

This simple analysis gives rise to a useful concept. Namely, the above flow of circles 
under the Riccati equation defines an invariant region consisting of the interior of the 
circles together with their boundary. In the next section we shall implement this point 
of view in the case of a non-constant potential V. 

3. Invariant Disk Estimates for the Riccati Equation 

We now turn our attention to solutions y G C of the Riccati equation (II. 4p for a 
general complex- valued potential V{x). Our goal is to derive inequalities of the form 

\y{x) — 7n{x)\ < R{x) 

for explicit functions m and R. Geometrically, this means that the exact solution y 
stays inside the disk of radius R centered at m. Thus we can interpret m as an 
approximate solution, with R giving a rigorous bound on the error. 

In preparation, we show that the the flow given by the Riccati equation maps circles 
into circles. At the same time, we derive differential equations describing the motion 
of these circles. Let us assume that y lies on a circle of radius R centered at m, i.e. 

y = m + Re'^P, 



ERROR ESTIMATES FOR THE RICCATI EQUATION 

Im y 





In 


y 




/^ 


7<~;> j 


I 




^ 




U\ 




• 




Vj 


^ 


-c 


0.5" 


■( 7~ 


\ 


-2 


-1 


; ^t/2 


/\ 




-0.5 


' ^5 


C/%. 




-1.0 


; 


"^Z 




» Rey 



Figure 1. Example for flow of circles with C, = 2 — i (left) and station- 
ary circles in the case C = ^ (right). 



where i? > 0, (/? E M, ttt, G C. We now derive differential equations for K and ra. First, 

m! ^ b! e^ ^ i^' ne^ = y' 

= V -y^ = V -m^ - 2Rme''f - R^e^'"^ . 

We multiply by e"**^ and take the real parts to get 

Re [e-''^{V -m^ - m') - 2Rm - R' - iip' R - e'"^ R^] = . 

In taking the real part, the term involving c^' drops out. Furthermore, we may replace 
the factor e**^ by e"**^. This gives 



Rele-'-^iV 



m 



m' - R^) - 2Rm - R'] 



0. 



Since this equation should hold for any (/?, we obtain the system of ODEs 

R' = -2RRem 

m' = V -m^ -R^ . 

It is difficult to solve this system exactly. But the following lemma gives us a method 
to obtain invariant disk estimates. 

Lemma 3.1. For a given complex-valued potential V = V{x), we assume that the 
complex function ni and the real function R > satisfy on the interval I := [xo,a;i] 
the differential equations 

Rl = -2RRem + 6R (3.1) 

m' = V -m'^ -R^ + 6m, (3.2) 

where the real function 5R and the complex function 5m obey the inequality 

5R > \6m\ . (3.3) 

Then the disk with radius R centered at m is invariant, i.e. for any solution of the 
Riccati equation (|1.4p the following implication holds: 

\y — 7n|(xo) < R{xo) =^ \y — m\{x) < R{x) Vx G / . 
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Proof. We first prove the lemma in the case of strict inequality in (j3.3p . If the lemma 
were false, there would be an x € / such that \y — m\{x) > R{x). Let x G [xo,xi) be 
the infimum of all such x. Then \y — m\{x) = R{x). Using (jl.4p and (|3.2p . we find at x 



d . , ^ (y' - m')(y - m) ^ (V - y"^ - V + m^ + R^ - 5rn){y - m) 
'y — m\ = Re j j = Re 



dx \y — m\ \y — m\ 

< fiRe(,-m) + |M-Re ^^'-,"^'^(^,-^ 

\y — m| 

= RRe{y — m) + \6m\ —RRe^y + m) = —2RRem + \6m\. 
Combining this with (13. ip and the strict inequality in (13. 3p . we obtain 

-— |y-m|(;r) < --^R{x) . 
ax dx 

Hence there is a 5 > such that on the interval (x, x + 5) the inequality \y — m\ < R 

holds. This is a contradiction. 

To prove the general case of non-strict inequality in (j3.3p . we introduce the new 

radius R by 

R = i? + ee^(^-^'i), 

where e > and 

L > max('-2Rem + 2i?-ee^(^-^i)V (3.4) 



Then R and m satisfy again the differential equations (13. ip and (13. 2|) . where R is 
replaced by R, and 6R and 6m are replaced by 

6R = (5i? + ee^(^-^i)(L + 2Rem) 

dm = 5m + ee^(^-^i)f2i?-ee^(^-^i)' 



Using (j3.3p and (|3.4p . we see that |5rfi| < 5i? for sufficiently small e, and thus the circle 
with radius R and center tti is invariant. Letting e \ completes the proof. D 

In what follows, it is convenient to write m = a + ifi with real a and /?, and similarly 
we write 5m = 5a + i 5(3. Then the equations ()3.ip and ()3.2|) become 

i?' = -2afi + 5/? (3.5) 

a' = ReV - a"^ + (3^ - R^ + 5a (3.6) 

(3' = ImV -2aP + 5/3. (3.7) 

In our first theorem we satisfy the conditions of Lemma 13.11 by employing a special 
ansatz where the invariant disks are centered on the real axis (so that (3 = 0). We 
introduce the abbreviation 



U = ReU-a^-a 



(3.8) 



Theorem 3.2. Suppose that on an interval I C M, we are given the real func- 
tions a,W (z C^{I). Assume furthermore that W > and that 

2aW + — - VW (\W - U\ + \ImV\j > 0. (3.9) 

Then the circle with radius R centered at m with 

R = ^/w , m = a (3.10) 
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is invariant on I under the Riccati flow (jl.4p . 

Proof. We shall verify the conditions of Lemma 13.11 Using the ansatz (|3.10p in the 
equations (jS.Sp . (j3.6p and (j3.7p . we get 

6R = -^ + 2aVw, 6a = W-U, 5(5 = -ImF. 

2^/W 



The sufficient condition 5R > \5a\ + |5/3| reduces to (13.9J) . D 



In our next theorem we will not prescribe (3. For convenience, we now replace (13. 3p 
by the sufficient condition 

(5i? = |5a| + \5f3\ . (3.11) 

Thus our task is to find six functions i?, a, /3, 5R^ 6a and 5/3 which satisfy the 
equations ()3.5p . (13. 6p and ()3.7p together with the condition ()3.1ip . Thus we expect 
that we can freely assign two functions, and then the remaining four functions will 
be determined. To explain which functions are to be preassigned, we note that the 
equations ()3.5p and (13. 7p are both linear, so that the only nonlinearity appears in the 
a-equation ()3.6p . Thus we can get rid of the nonlinear differential equation by taking 
the functions 

a and W := R^ - (3'^ (3.12) 

as a-priori given. Then (j3.6p becomes the defining equation for 6a] namely, 

6a = -Rey + a^ + a' + W . 
Using these definitions, we are left with the system of equations 

R! = -2aR + \6a\ + \6l3\ (3.13) 

p' = lmV-2af3 + 6p (3.14) 

for the three unknown functions R, (5 and 6(3. This system consists of two linear 
differential equations and one nonlinear algebraic equation. Of course, we must also 
ensure that < i? < cxd. 

To analyze these equations, we consider the two cases 5/? > and 5/3 < separately. 
Case (A) 6[5 > 0: Subtracting (fXTip from (I3l3]l we get 

{R-py = -2a{R-/3)-lmV+\6a\. (3.16) 

Furthermore, differentiating (I3.15P gives 

2RR'-2PP' = W' . (3.17) 

We consider (|3.16p and (j3.17p as a system of two linear algebraic equations in 
the unknowns R' and (3'. Solving for (3' and substituting into (j3.14p . we find 

W 
(R -(3)6p = -R \6a\ + 2Wa + ImV (^ + — . (3.18) 

Since we assumed 6(3 to be positive, we get the consistency condition 



{R-(3) 



W 
2aW + — -R \6a\ + [3 ImV 



> 0. (3.19) 



Provided that this consistency condition holds, we can use (j3.16p to solve 
for R- p. Then using (|3T5]) . we conclude that R + f3 = W/{R - (3). 
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Case (B) 6(3 < 0: Adding iKUh and (|XTij) gives 

{R + py = -2a {R + (3) + lmV + \6a\. (3.20) 

Solving (j3.20p and (j3.17p for /?' and substituting the result into (j3.14p . we 
obtain 

w 

{R + (5)5I3 = R \5a\ - 2Wa - IxnV P-^ , 
giving rise in this case to the consistency condition 



{R + (3) 



W 
2aW + — -R \5a\ + [3 ImV 



> 0. (3.21) 



Provided that this consistency condition holds, R + (3 is determined by (j3.2Up . 
and R-P = W/{R + P). 

Since the bracket in (j3.19p and (j3.2ip will appear often, it is convenient to introduce 
the abbreviation 



W 
D = 2aW + R \6a\ + (3 ImV 



(3.22) 



We refer to 2) as the determinator. Then the following invariant disk estimate holds. 

Theorem 3.3. Suppose that on an interval I C M, we are given the real func- 
tions a,W £ C^{I). We set 

a{x) = exp ( / 2a ) . (3.23) 



Then the following statements hold. 

(A) Define real functions R and P on I by 

{R-I3){x) = - I (T{-lmV + \W -U\) (3.24) 

{R + P){x) = ,J^^':l , (3.25) 

(R- P)[x) 

(with U as defined by ()3.8p ). Assume that the function R — (3 has no zeros, 
R> 0, and that 

{R-P)'D > (3.26) 

(with 2) according to (I3.22p ). Then the circle centered at m{x) = a + i[3 with 
radius R{x) is invariant on I under the Riccati flow (|1.4p . 

(B) Define real functions R and (3 on I by 

{R + P){x) = - f a {ImV + \W - U\) (3.27) 

Assume that the function R + (3 has no zeros, R> 0, and that 

{R + P)Ti > 0. (3.29) 

Then the circle centered at m{x) = a + i(3 with radius R{x) is invariant on I 
under the Riccati flow ()1.4p . 
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Proof. We only prove the first part, as the proof of the second part is similar. Since (j3.26p 
holds, we are consistently in Case (A) above. Integrating the differential equation (j3.16p 
gives (j3.24p , whereas (|3.25p follows from (j3.12p . We thus have a solution of the system 
of equations (|3.13p , (|3.14p and (|3.15p with 6(3 given by (|3.18p . Applying Lemma 13.11 
completes the proof. D 

We next outline the general strategy for applying the last theorem. Suppose that 
an approximate solution y of the Riccati equation is given. Setting V = y' -\- y'^ , we 
can consider y as a solution of a Riccati equation 

y' = V-f (3.30) 

with an "approximate" potential V. Since y approximates y, an obvious idea would 
be to consider circles centered at m = y. In the estimates of Theorem 13. 3|, however, 
we only need to prescribe the function a = Rem, whereas the imaginary part of m 
will then be determined. Therefore, it is natural to choose 

a = Rey . (3.31) 

A natural choice for W is to set 

W = U (3.32) 

(with U according to (j3.8p ). because then the term |T^ — C/| in (j3.24p and (j3.27p 
vanisheqj. Then Theorem 13.31 gives us invariant disk estimates. The center m of this 
disk can be viewed as an improved approximate solution, and its error, measured as 
the distance from the exact solution, is at most R. The invariant disk estimate depends 
crucially on the sign of the determinator 1). Thus in applications it is important to 
know the sign of 2). It can even be useful to modify y so as to give 2) a desired sign 
(this will be illustrated in Section [5] in a few examples). For this purpose, it is useful 
to bring the determinator into a more convenient form. 

Lemma 3.4. Suppose that, for a given solution y of an approximate Riccati equa- 
tion (I3.30p . a and W are chosen according to (j3.3ip and (j3.32p . Then the determinator 
is given by 

D = 2aRe{V-V) + -Re{V -V)' - pimV + pimV , (3.33) 

where (3 = Imy. 

Proof. Writing the Riccati equation (|3.30p as 

a' = Rey-a2 + ^2^ ^' = Imy-2a/3, 

we can use p.32p and (|3.8p to obtain 

W = U = Re(y - y) - /3^ 
Applying these identities in (j3.22p gives 

S) = 2aR(i{V -V) + -R(i{V -V)' -2a^'^ - ^^' + (3luiV . 

Using the /^'-equation finishes the proof. D 



Other choices of W might be useful, but we will not consider them in this paper. 
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4. Error Estimates for Real Potentials and Examples 

In this section we shall illustrate our error estimates in examples for a real poten- 
tial V and explain how the previous results from [3l [2] fit into our framework. We 
first mention an additional structural relation which appears only for a real potential. 
Namely, suppose that (/) is a complex solution of the Schrodinger equation (jl.ip . De- 
composing it into its real and imaginary parts, (j) = (j)i+i(p2, the functions c^i and (j)2 are 
real solutions of the Schrodinger equation forming a fundamental set. Thus their Wron- 
skian w = 4>i4>2 — 4'i4'2 is a constant. It can be expressed as w = Im ((/)(/>') = |</>p Imy, 
giving rise to the relation 

Im y 
which is useful for estimating the amplitude of (p. Moreover, this relation shows that 
the imaginary part of y can never change sign, and thus the upper and lower half 
planes must be invariant regions for y. 

To illustrate how our estimates apply, we begin with a simple example. 

Example 4.1. (negative, increasing potential) Suppose that the potential satis- 
fies on the interval / = [xo,Xi] the conditions 

V{x) < 0, V'{x) > 0. 

The simplest (but certainly not the best) choice is to set a = 0. Furthermore, set- 
ting W = U, we see from (|3.23p that cr is a constant. According to (|3.8p . we find 
that W = V. The determinator (j3.22p simplifies to 

V 
2) = - > 0. (4.1) 

We want to apply part (B) of Theorem 13.31 (part (A) is similar if one flips the sign 
of 13). From (I3.27P we see that R + (3 is a constant, which according to (I3.29P and (14. ip 
must be positive. We thus obtain 

p + R = c>0, (3-R=^-^. (4.2) 

c 

The condition R> gives the constraint 

c^ > max|y| = \V{xo)\ . 

If the potential is constant, we recover the stationary circles of Section [2] (as is easily 
verified by comparing (|4.2p and (|2.ip l. If the potential is non-constant, we get an 
increasing family of invariant disks in the upper half plane centered on the imaginary 
axis whose highest point is fixed at ic, and whose lowest point varies like i|y|/c; see 
Figure El 

We next consider the situation for a general real potential V and a general choice 
of a, again setting W = U. Then the determinator becomes 

ir_ _ u{aHjy 

T ~ 2" 0-2 f/ 

(with a defined by (|3.23p ). The main simplification is that now £) is an a-priori given 
function. In the case when R — (3 and R + (3 have opposite signs, we see from ()3.26p 
and (j3.29p that, depending on the sign of 2), either part (A) or part (B) of Theorem l3.3l 
applies. In this case, we can combine parts (A) and (B) to obtain a statement which 
can be used even if 2) changes sign infinitely often. In view of (j3.15p . the condition 



S = 2af/ + - = -^-^rf (4.3) 
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Figure 2. A simple invariant disk estimate for a negative, increasing potential. 



that R — (3 and R + (i should have opposite signs is equivalent to the simpler condition 
that [/ < 0; it means that the invariant disk either lies in the upper half plane or in the 
lower half plane, without intersecting the real line. By continuity, the invariant disk 
cannot move from the upper to the lower half plane. Thus we can assume without 
loss of generality that the invariant disk lies in the upper half plane. This means 
that R + (3 > Q and R — (5 < 0. It is convenient to satisfy the equation (|3.15p (and 
thus also (13.251) and p.28p ) by making the ansatz 



P 



T.i 



R 



T 



2 \ TJ 2 

with a free function T > 1. If 2!) is negative, we apply Theorem 13.31 (A) to obtain 



(4.4) 



R-P c 

with c an integration constant. Diffentiating gives 



r 



i'^vm' 



a 



1 {a^\u\y 

2 



1 {a^uy 

2^W 



{a^uy 



a^U 



'\U\ 2 a^\U\ 

where in the last two steps we used that U is negative, and that in the considered 
case 2) < 0, the term {a'^Uy /{a'^U) is positive according to ()4.3p . If 2) is positive, we 
apply Theorem 13. 31 (B) to get 

R + P c 



Now differentiating gives 

r _ _(ct 



\u\y 



a 



\U\ ay/\U\ 

1 {a^uy _ 

2 a^U 



{a^uy 



a'U 



where in the last step we used that now {a'^Uy / {a'^U) < 0. We conclude that, inde- 
pendent of the sign of 2), we obtain the differential equation 
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Figure 3. Invariant disk estimate for U < 0. 



Integrating both sides, we can rewrite T as a total variation, 



T{x) = To exp Q TV [,„,,) log \a^U\] 



(4.5) 



In this way, we recover the invariant disk estimate which was first obtained in [SJ 
Lemma 4.1]; see Figure [3l 

Lemma 4.2. Let a be a real function on [xq,xi] which is continuous and piecewise 
C^ . Suppose that the function U := V — a^ — a' is negative on [xq^xi]. Using the above 
definitions (]4.5p and (j4.4p , the disks of radius R centered at m = a + i/3 are invariant 
under the Riccati flow (II. 4[) . on [xo,xi]. 

We now illustrate in two examples how this lemma can be applied (for other ap- 
plications see 1^ Theorems 4.2 and 4.3] and [2, Lemmas 4.10 and 4.12]). In the first 
example, we consider the case that the potential V is negative. We recall that the 
WKB wave functions are given by (see for example [21 Section 2.4] or [5]) 



0wkb(2;) 



\V{x)\~-^ exp (±i I y/\V\) ■ 



(4.6) 



A short calculation shows that these functions satisfy the equation 

dV ... /5 |y|'2 1 \V\"^ 



^-^(AwKB - V<t)- 



WKB 



16 \V\^ 



4 \V\ 



^WKB 



Interpreting the right side as an error term, one sees that the WKB wave functions 
are expected to be a good approximation to (j) provided that 

y'2 



V" 



V 



y3 



« 1. 



(4.7) 



For definiteness, we consider the plus sign in (j4.6p (the minus sign can be obtained by 
complex conjugation). Then the corresponding approximate Riccati solution is given 

by 

V 

2/WKB 



^WKB 



W 
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Following the general strategy described after Theorem 13.31 we thus choose 



v 

a = ReywKB = -777 ■ (4. 



An easy computation gives 



V" V^ 

/ V" 5y'2 \ 



In view of (j4.10p . the condition (j4.7p and the fact that 1/ < 0, we see that U is negative. 
Thus Lemma 14.21 applies, giving the following estimate. 

Example 4.3. (error estimate for WKB wave function, y < 0) 

Suppose that V is negative and 

~4y2 ~ I6|y|3 < ^• 
Then the disks of radius R centered at m = a + ij3 as given by ()4.8p . ()4.4p and 



T{x) = To exp - TV^xo.x) log 



1 



2 



V" sy 



/2 



1 + TT7T7 + 



41/2 i6|y|3 



(4.13) 



are invariant under the Riccati flow (|1.4p . on [xo,xi] 



In our setting, the semiclassical limit can be described by scaling the potential 
according to y ^ A^ and considering the behavior as A ^> 00. In this limit, the square 
bracket in ()4.13p converges to one, so that the logarithm vanishes. As a consequence, 
T becomes a constant. Choosing Tq = 1, we find that the radius R{x) as given by (14.41) 
vanishes identically in the semiclassical limit, in agreement with the fact that the WKB 
wave function goes over to the exact solution. 

The next estimate even applies near a zero of the potential. 

Example 4.4. (invariant disk estimate with exponential bound) 

Choose a constant a, 

a = c+ sup v^max(0, V{x)) (4.14) 

x£[xo,xi] 



with c > 0. Then 



U = V-a^ < -c^ < 



a^U = C'' e^'^^ {y - a") 

with an integration constant C. Hence Lemma \4-^ yields that the disks of radius R 
centered at m = a + i(3 as given by ()4.14p . ()4.4p and 

T{x) = n exp Q Ty[,„,,) log [e^"- [V - a')]^ 
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Rey 



a + VU 



Figure 4. Invariant lens-shaped region estimate for [/ > 0. 



are invariant under the Riccati flow (II. 4p . on [xojXi]. 

Clearly, due to the exponential factor e^°^, this estimate gets weaker as the size of 
the interval increases. We shall below see a much better estimate using Airy functions 
(see Example 15. ip . 

We next consider how Theorem 13.31 applies in the case f7 > (where we again 
choose W = U). Then, according to (j3.15p . the functions R — (3 and R + 13 have the 
same sign, which must be positive because i? > 0. Inspection of ()3.26p and (j3.29p 
yields that the determinator must be positive, i.e. 



U' + AaU > 0. 



(4.15) 



Provided that this condition holds, we may apply both parts (A) and (B) of Theo- 
rem [331 Taking the intersection of the corresponding invariant disks gives the following 
result, which generalizes [21 Lemma 4.2]. 

Lemma 4.5. Let a he a real function on [xq,xi\ which is continuous and piecewise 
C^ . Suppose that the function U := V — a^ — a' is positive on [xq,xi\ and that the 



condition (I4.15P holds. For any positive constants ct, k 
disks of radii Rk centered at ruk = a + i(3k with 



1,2, we introduce the two 



i?i 



R2 



Ua 
ci 

C2 



a 



/5i = ^ — -- 



P2 



1 fUa 

2 Vci 
1 f Ua 

'2\^ 



C2 



Then the two disks as well as the lens-shaped region obtained as their intersection are 
invariant under the Riccati flow (jl.4p . 

As is easily verified, the boundaries of both disks intersect the real axis at the points 
a lb VU, and thus their intersection is non-void; see Figure JH 
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In our next example we shall apply Lemma 14.51 to the WKB wave function for 
a positive potential. We thus assume that y > and that (j4.7p is again satisfied. 
Considering the WKB wave function 

</'wkb(x) = V{x)--^ expff Vv\ , (4.16) 

an easy computation shows that 



^WKB 



ywKB(x) = x^ = Vy-^ 

<PWKB 41/ 

r- y 

a = ReywKB{x) = V F - — 

U = v-a^-a' = --^ + — . 

16 V^ AV 

Unfortunately, in contrast to our previous example (j4.10p , now U need not have a fixed 
sign, and thus one cannot apply the above lemmas. Our method for getting around 
this problem is to consider in the WKB ansatz another potential Vwkbj and to set 

4>wkb{x) = VwKB(a;)"4 exp ( / \AVkbJ • 



Then 



ii 



K, 



a 



Rep™ = T^Vkb-:^^?^ (4.17) 

<Pwkb 4 Kwkb 

U = v-a'-a' = V-V^^^-^^^ + ^^. (4.18) 

16 KvKB 4V/wKB 

By arranging that U has a definite sign, we may apply either Lemma [4.2l or Lemma [4. 5 [ 
To give a simple example, we choose 

V 

VwKB — -r ■ 

Then U is positive because V is positive and (j4.7p holds. Furthermore, a direct com- 
putation gives 

a = exp( 2a\ = -^ exp ( / Wj (4.19) 



2 .. 3 

F^ = 2 + 

^ {U' + 4aU) = 1 + 



5 y'2 y" 
+ 



3yi 



8 y3 2^2 

^y^ 5^3 V^ WV" V" 



12 F3 ' 8 ^1 ' 3^2 ^^1 g^l 
We are now in a position to apply Theorem 13. 3[ 



(4.20) 
(4.21) 



Example 4.6. (error estimate for WKB wave function, y > 0) 

Suppose that V is positive and that the square brackets in (j4.20p and (j4.2ip are both 
greater than minus one. Then choosing U and a according to (14.201) and (14.19p . 
the lens-shaped region from Lemma \4-5\ (see Figure ^ is invariant under the Riccati 
flow ()1.4p . on [xo,xi]. 
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It is worth noting that from (j4.19|) . we see that a grows exponentiahy. Hence the 
lens-shaped region gets thinner, improving the estimate exponentially fast. 

The previous estimates using Lemma [4.2l and Lemma [4. 5 1 have the disadvantage that 
if \J changes sign, different kinds of estimates must be pasted together. This makes it 
necessary to match invariant disks (see [21 Lemma 4.4]). In the next section we explain 
a more convenient method in the more general context of a complex potential (but 
also working for real potentials) which makes it possible for the invariant disks to flow 
continuously across the real line. 

5. Error Estimates for Approximate WKB/Airy Solutions for Real or 

Complex Potentials 

In this section we illustrate the general technique for applying Theorem 13.21 and 
Theorem 13.31 bv discussing typical examples. Based on Theorems 13.21 and 13.31 we will 
derive rigorous error estimates for the standard approximate solutions obtained by 
glueing together WKB and Airy wave functions. 

For a given potential (real or complex), one can usually distinguish regions where \V\ 
is large, so that the solutions of the Schrodinger equation are well-approximated by 
WKB wave functions. In the remaining regions, one can approximate the potential by 
a linear potential, so that the Schrodinger equation can be solved explicitly in terms of 
Airy functions. To be more flexible, in the WKB region we consider the approximate 
solutions 

(?!>wkb(2;) = Vwkb(2;)~4 exp f / ±v^VVkb j (5.1) 

with a potential Vwkb ~ y, which will be determined later. Likewise, in the Airy 
region, our approximate wave functions are solutions of the Schrodinger equation 

<Pl{x) = Va<Pa, (5.2) 

where Va ~ T^ is a linear function, also to be determined later. By the standard C^- 
glueing of the WKB and Airy functions, we obtain an approximate wave function 0. 
The wave function is a weak solution of a Schrodinger equation 

4>" = V4>, (5.3) 

where V in the Airy region coincides with Va, whereas in the WKB region a short 
calculation gives 

V — V -\- WKB '-^WKB (c- A\ 

16 V^KB 4FwKB 

Note that V is piecewise smooth, but in general has discontinuities. Following the 
strategy explained after Theorem 13.31 we introduce the function y = </>'/0, which is 
continuous and piecewise smooth, satisfying the Riccati equation 

y' = V-f . (5.5) 

We again introduce the (continuous, but only piecewise smooth) function a by 

a = Rey , (5-6) 

define U by (|3.8p and (again for simplicity) choose W = U. We choose starting 
values for the functions (3 and R which are compatible with (I3.15[) . We can then 
apply Theorem 13.31 in various ways. If C/ < 0, i.e. R + (3 and R — (3 have opposite 
signs, depending on the sign of the determinator 2), we can apply either part (A) 
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or part (B). If however U > 0, the functions R + j3 and R — (5 are both positive, 
and thus Theorem 13.31 apphes only if 5} is positive, and in this case we can apply 
both part (A) and part (B). If (3 becomes zero, it may be preferable to switch to the 
simpler estimate of Theorem l3.2| provided that the inequality (j3.9p holds. In this way, 
we have different possibilities to obtain invariant disk estimates, and by taking their 
intersection one gets even sharper estimates involving lens-shaped invariant regions. 
Observe that on the boundaries between the WKB and Airy regions, the functions V 
and a', and consequently also [/, may have discontinuities. In order to satisfy (j3.15p . 
either R + (3 or R — (3 must "jump" discontinuously, in such a way that the new disk 
contains the old one. 

The crucial point for making this procedure work is that one must be able to pre- 
scribe the sign of 2) or to satisfy the inequality (j3.9p . One method to achieve this is 
to modify the potentials Vwkb and V\ in (j5.ip and (j5.2p . The effect on 5} of modify- 
ing the potentials can be seen most easily from Lemma 13.41 More specifically, in the 
WKB region, the term 2aRe(y — V) in ()3.33p can be suitably changed by modifying 
the real part of Vwkb; whereas the imaginary part of Vwkb affects the term —(51vi\V. 
Note that the derivative term in (J3.33P as well as the derivative terms in (15.41) only 
give rise to small corrections. In the Airy region, on the other hand, we can change 
the coefficients in the linear function Vx{x) to modify the suitable terms in (j3.33p . In 
order to satisfy the condition (13. 9p . it may be useful to apply the identity 



dSj]) = 2aRe{V-V) + ]- Re{V -V)' -^luiV - \lmV\\jRe{V - V) - (3^ 

(keeping in mind that in this case, /3 = 0, and thus f3 should be small), or else one can 
modify a, without respecting the relation ()5.6p . 

We now give two concrete examples to illustrate this method. In the first example, 
the imaginary part of the potential is so small that the method works just as well for 
a real potential. 

Example 5.1. We consider the potential (inspired by the spheroidal wave operator [3j) 

V = 10000 (-- + (1-F 0.05 i) sin^x j , (5.7) 

on the interval [0, ^]. Due to the large prefactor, the WKB condition (j4.7p is satisfied 
except in a small neighborhood of the point x = j, where \V\ is small. Thus we divide 
our interval into three regions: 

(a) The "classically ahowed" WKB region [0,0.715] 

(b) The Airy region [0.715, 0.83] near the "classical turning point" 

(c) The "classically forbidden" WKB region [0.83, f ] 

(The precise choice of the boundary points is arbitrary and has no major effect. More- 
over, since the factor 0.05 is so small, we simply adopted the terminology from quantum 
mechanics, disregarding the effect of the imaginary part of the potential.) In the Airy 
region we approximate the potential by its linear Taylor series around the point x = j. 
We glue together the corresponding WKB and Airy functions, starting at x = with 
the WKB wave function (15. 1[) with Vwkb = V and choosing the plus sign. In the 
region (c) we also choose Vwkb = V- This gives the approximate wave function (j), 
and we define the corresponding approximate Riccati solution by y = (j)' /(j). Following 
the general procedure described earlier in this section, we choose a according to (j5.6p 
and set W = U. To better illustrate our estimates, we choose the radius R of the 
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Figure 5. Upper and lower bounds for the imaginary (left) and real 
(right) parts of the Riccati solution in Example 15. 1[ 



initial disk to be 2.5, although choosing it equal to zero would give a better estimate. 
In Figure [5] (left), the imaginary part of an exact numerical solution y starting on 
the boundary of the initial circle is plotted together with the upper bound (3 + R and 
the lower bound (5 — R. Likewise, in Figure O (right), the real part of the solution as 
well as the upper bound a + R and the lower bound a — R are given. On sees that 
the invariant disks jump discontinuously at the glueing points, such that the new disk 
contains the old disk. In Figured! 

the invariant disks are plotted for discrete values of x, and the black dots denote 
the exact solution at these values of x. We point out that the invariant disk always 
stays in the upper half plane. In region (c), the determinator is negative, so that part 
(A) of Theorem 13.31 applies. This means geometrically that the lower bound [i — R 
approaches the exact solution exponentially fast, whereas the upper bound (5 + R \s 
not as good an approximation. 

In order to get a better upper bound, we changed the WKB potential in region (c) 
to 

^WKB = 0.9 y. (5.8) 

This makes the determinator negative, so that part (B) of Theorem 13.31 applies. 
Since Vwkb in (jS.Sp deviates considerably from V , we cannot expect that our ap- 
proximation will be good in region (c). This is reflected in Figure [7] by the fact that R 
becomes large if x approaches ^. Nevertheless, this estimate is still useful because the 
upper bound P + R approaches the exact solution exponentially fast. Taking the inter- 
section of the invariant disks with those from Figure O one gets lens-shaped invariant 
regions, thus estimating the exact solution up to exponentially decaying errors (see 
Figure E]). This exponential shrinkage of the invariant region can be understood by 
the fact that the the potential is slowly varying, and that for a constant potential the 
solutions would tend exponentially fast to the stable fixed point \/V (see Section [2]). 

Finally, in Figure [9] we show invariant disk estimates, where in region (c) we applied 
Theorem 13.21 and chose f3 = 0. This estimate is not as good as the previous estimate, 
but on the other hand Theorem 13.21 is easier to apply, and the coarser estimate might 
be sufficient for some applications. 

In the previous example, the semi-classical approximation was very good because 
of the large prefactor in (j5.7p . In our next example, the WKB/Airy wave function is 
not such a good approximation, making the error estimates more subtle. Moreover, 
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Figure 6. Invariant disks for the Riccati solution in Example 15. 1[ 
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Figure 7. Upper and lower bounds after flipping the sign of S) in Example 15. 1[ 




Figure 8. Lens-shaped region estimates for the Riccati equation in Example 15. 1[ 



the imaginary part of V has the opposite sign, thereby forcing the Riccati solution to 
cross the real axis. By suitably adjusting the sign of the determinator, we will arrange 
that the invariant disks also move from one half plane to the other. 
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Figure 9. Invariant disk estimate applying Theorem 13.21 in Example 15.1 



Example 5.2. We consider on the interval [0, ^] the potential 



V 



500 



+ (l-0.2i)sin2x 



We choose the WKB regions to be [0,0.52] (region (a)) and [0.83, |] (region (c)), 
and the intermediate region [0.52,0.83] (region (b)) is the Airy region. In Figures [TU] 
and [TT] upper and lower bounds are given. In Figure [TJ] we plot the corresponding 
invariant disk estimates as well as lens-shaped invariant regions obtained by taking 
the intersection of the disks of Figures [10] and [TTl 

Note that now the invariant disks cross the real axis. When this happens, one must 
ensure when applying Theorem 13.31 that the denominator in (I3.25P or (I3.28P never 
becomes zero, because otherwise we would lose control of the estimates. This means 
that if the function R — (3 vanishes, we must arrange the sign of the determinator such 
that case (A) applies, whereas R + f3 may vanish only when we are in case (B). To this 
end, the potentials Vwkb and Va are chosen as follows: 




Figure 10. Upper and lower bounds for the imaginary (left) and real 
(right) parts of the Riccati solution in Example 15. 2[ 
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Figure 11. Upper and lower bounds after flipping the sign of 2) in Example 
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Figure 12. Invariant disk (red) and lens-shaped region estimates 
(blue) for the Riccati equation in Example [ 



In region (a), we choose Vwkb = V. In region (b), for Va we took the linear Taylor 
polynomial and decreased the real part of the linear term so as to make 1) positive. 
Then we are in case (B) of Theorem l3.3| in which the function f5 — R can cross smoothly 
across the real line. Shortly after [5 — R has become negative, we switch to case (A) of 
Theorem 13.31 so that (3 + R can smoothly flip sign. For the estimates of Figure [10] we 
choose the potential Vwkb = V^ in region (c) such that 2) is always positive. Then f3—R 
is a good approximation. For the estimates of Figure [TTl however, we modified Vwkb 
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Figure 13. The potential V — V as used for the estimates in Figure [TO] 
(left) and Figure [H] (right) . 



in region (c) such as to arrange that S becomes negative, with the result that (3 + R 
is a good approximation. 

Our choice of the approximate potentials is shown in Figure [T3t where we plot the 
real and imaginary parts oiV — V , with V as in ()5.3p and (15. 5p . We point out that the 
scale in these plots is one order of magnitude smaller than the scale of y ~ 250. This 
illustrates that slight modifications of the approximate potentials suffice to arrange the 
appropriate sign of 2) . The detailed form of the choice of the approximate potentials 
affects the function /3 + i? in Figure [10] (left) or the function /3 — i? in Figure [TT] (left), 
as well as the functions a it i? on the right of these figures. However, the corresponding 
lens-shaped invariant regions are insensitive to the detailed choice of V] see Figure [T2l 
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